Feasibility of a synchrotron storage ring for neutral polar molecules 
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Using calculations and mathematical modeling, we demonstrate the feasibility of constructing a synchrotron 
storage ring for neutral polar molecules. The lattice is a racetrack type 3.6 m in circumference consisting of 
two of 180-degree arcs, six bunchers, and two long straight sections. Each straight section contains two triplet 
focusing lenses and space for beam injection and experiments. The design also includes a matched injector and 
a linear decelerator. Up to 60 bunches can be loaded and simultaneously stored in the ring. The molecules are 
injected at 90 m/s but the velocity of the circulating beam can be decelerated to 60 m/s after injection. 

The modeling uses deuterated ammonia ( I4 N 2 H 3 ) molecules in a weak-field seeking state. Beam that survives 
400 turns (15 s), has horizontal and vertical acceptances of 35 mm-mr and 70 mm-mr respectively, and an energy 
acceptance of ±2%. 
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PACS numbers: 29.20 Dh, 41.75.Lx, 33.80.Ps, 39.90.+d, 33.55.Be 



I. INTRODUCTION 

An electric field gradient exerts a force on the dipole mo- 
ment of a neutral polar molecule. The force, F x , in the (trans- 
verse) x direction is: 



dW _ dWdE 
dx dE dx 



(1) 



where W is the potential energy, in an electric field of mag- 
nitude E = (Ex+Ey) 1 / 2 (Stark effect), of the molecule. A 
similar expression may be written for the force in the y direc- 
tion. 

It is therefore possible to accelerate (decelerate) [|l| |[ |3]|. 
deflect, focus, and store polar molecules. Prior to this study, 
a continuous torroidal ring, without a buncher or matched 
injector was suggested by Auerbach et al.[Q], analyzed by 
Katz[Qj, and Loesch and Schell [^|, and recently constructed 
by Crompvoets et al.[[7]], who captured single pulses (about 
10 6 molecules) of deuterated ammonia ( 14 N 2 H 3 ) at 90 m/s (ki- 
netic energy w 9.7 K) and observed them for six turns (0.053 
s). 

A ring with much longer storage times that can capture and 
store the high peak intensity of (the decelerated output of) a 
pulsed molecular beam jet source, could increase the avail- 
able flux. With bunching, the energy of the stored beam can 
be varied continuously and the density of the molecules var- 
ied. These features make a molecular synchrotron storage ring 
useful for high-resolution spectroscopy, low energy scattering 
experiments, and for evaporative cooling [|| S], With evapo- 
rative cooling, the molecules can reach ultra-low temperatures 
where they may form quantum condensates. 



*E-mail H"Nishimurii@lhl.gov 
^E-mail: GRLambertson @ 1 b 1 .gov 
*E-mail: IGKalninsS'lbl.go' 
§ E-mail: HAGould@lbl.gov 



Evaporative cooling requires: high densities to thermalize 
the molecules by elastic scattering, forced evaporation of the 
hottest molecules, and long storage times for many repetitions 
of the cooling cycle. Bunching (and focusing) raises the den- 
sity to thermalize the molecules and beam drift (and/or de- 
bunching) isolates the hottest molecules for removal. 

In this paper we present the results of a study showing 
the feasibility of constructing a molecular synchrotron storage 
ring that will capture large numbers of molecules in multiple 
bunches, maintain the bunches, store the molecules for longer 
than 15 seconds and vary their kinetic energy. 

The paper is organized as follows: In section |l| we dis- 
cuss the lattice, the choice of molecular state, the necessity of 
avoiding regions of very weak electric field, and the dynam- 
ics calculations. In section |U]| we describe the performance 
of the synchrotron storage ring, including the dynamic aper- 
ture, operating parameters, the effect of gravity, bunching and 



deceleration, and collisional losses. Finally, in section IV 



evaluate the overall performance of the storage ring by adding 
a molecular beam source, a linear decelerator, and an injection 
system. 



II. STORAGE RING DESIGN 
A. Lattice 

Many different designs and placements of the elements 
were considered and analyzed. Our resulting design (Fig.[j]) 
is a racetrack configuration consisting of two of 180-degree 
arcs, 4 triplets of focusing lenses, six bunchers placed in re- 
gions having low horizontal dispersion, an injector for loading 
up to 60 bunches into the ring, space in the long straight sec- 
tions for an extractor (for example a mirror image of the injec- 
tor), and for experiments or evaporative cooling. The lengths 
of the elements and drift distances between them are listed in 
Table I. The strong combined horizontal and vertical focus- 
ing in the bend sections makes the ring very compact. The 



FIG. 1: Schematic diagram of the synchrotron storage ring mod- 
eled in our study. The bend radius is 0.2 m and the circumference is 
3.36 m. QF and QBF are horizontal focusing (vertical defocusing) 
lenses, QD are horizontal defocusing (vertical focusing) lenses and 
BU are bunchers. The injector and a portion of the decelerator are 
also shown. Additional details of the lattice are found in Table I, de- 
tails of the focusing lenses and bend elements are found in Table II 
and Fig. 's | 5| ^ . D etails of the decelerator and injector are found in 
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TABLE I: Element length and placement for one-fourth of the lattice 



Element 



Length (cm) 



Cumulative travel (cm) 



Drift 


16.0 


16.0 


Buncher (BU) 


2.0 


18.0 


Drift 


2.0 


20.0 


QF 


6.0 


26.0 


Drift 


3.0 


29.0 


QD 


8.0 


37.0 


Drift 


3.0 


40.0 


QFB 


7.0 


47.0 


Drift 


4.0 


51.0 


90° Bend 


31.4 


82.4 


Drift 


0.5 


82.9 


Half buncher (BU) 


1.0 


83.9 



triplet of straight-lens electrodes then transforms the beam to 
a wider, more collimated beam that can traverse the 32-cm- 
long straight section as shown in Fig. |(a). The large varia- 
tion of focusing strength is apparent in Fig. ||(b), the plot of 
the focusing parameters (3 X and p v around the ring. 

The horizontal betatron phase advance in the 180-degree 
bending region has been made equal to 4k in order to have 
zero dispersion at the bunching electrode, as shown in Fig. ^| 
The very low dispersion at the bunching electrodes prevents 
the development of strong synchro-betatron coupling and pre- 
serves the momentum acceptance. 



B. Molecular state 

Many different molecules and kinetic energies could be 
used for our feasibility study. We selected 90 m/s 14 N 2 H3 
molecules in the weak-field seeking upper-inversion level 





distance (m) 



FIG. 2: Beam envelope (a) and focusing parameter P (b) in the hor- 
izontal (x) and vertical (y) directions. (3 is the distance in which the 
transverse (betatron) oscillation advances in phase by one radian. A 
schematic of the lattice is shown above for location reference. 
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FIG. 3: Horizontal dispersion of the beam around the lattice. Low 
dispersion at the bunchers prevents the development of synchro- 
betatron oscillations. 



(lowest vibrational state) \J,K,M >= [1,1, — 1 > (Fig. f|) be- 
cause they had been previously used by Crompvoets et al.|Q], 
and because the focusing properties of weak-field seeking 
states make the design and feasability assessment easier and 
the ring more compact. 

Molecules in weak-field seeking states (molecules whose 
potential energy increases in an electric field) can be fo- 
cused simultaneously in both transverse directions, whereas 
molecules in strong-field seeking states can be focused in only 
one transverse direction while defocusing in the other trans- 
verse direction. 

Fringe fields seen by the molecules upon entering and exit- 
ing the bending, focusing, and bunching elements will also ex- 
ert transverse forces on the molecules. A molecule in a weak- 
field seeking state traveling in the z direction will experience 
a net focusing effect in the vertical {y) direction upon entering 
or exiting an electric field produced by horizontal (x,z) plane 
parallel electrodes. A molecule in a strong-field seeking state 
will defocus. 

Except for the rotational state 7 = (which is always strong 
field seeking), each rotational state, J, contains M components 
which are partially or completely degenerate in zero electric 




2 4 6 

electric field MV/m 

FIG. 4: Stark effect in the \J,K >= |1, 1 > levels of 14 N 2 H 3 . The 
M = — 1 level is weak-field seeking and used for this study. Inver- 
sion splitting, dipole moment and rotational constant are taken from 
Townes and Schawlow [^o|.Hyperfine splitting and the Stark effect 
in the lfr 24 Joule inversion-split M = levels are neglected. 
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FIG. 5: Maps of (truncated) equipotentials in x and y of (a) F lenses, 
QF and QFB (horizontal focusing, vertical defocusing), (b) D lenses, 
QD (horizontal defocusing, vertical focusing), and (c) bend elements 
(horizontal and vertical focusing). The ellipse and circles show the 
approximate size of the beam envelope in the element. Note that the 
electric field direction at the orbit of the molecules is always in the 
vertical direction. The actual electrode can be fabricated to lie along 
any set of equipotentials that are larger than the beam size. 



TABLE II: Bending and focusing/defocusing electrodes 



Electrode 


Length (cm) 


E Q (MV/m) 


Half-gap (mm) 


A 2 


A 3 


bend 


31.4 


3.37 


6 


157 


5667 


QF 


6.0 


2.88 


15 





2000 


QD 


8.0 


3.55 


15 





-2000 


QFB 


7.0 


4.30 


15 
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FIG. 6: Electric field magnitude around the ring. All electric fields 
are in the vertical direction and do not change sign. The buncher field 
is shown fo r an o n-energy molecule entering or exiting the buncher 
(see section pic| ). The frin ge fiel ds are part of the focusing system 
and are discussed in section II D 3 



gradient sequence [|14||. To prevent rapid changes in elec- 
tric field direction, the field direction, at the molecule's orbit 
throughout the ring, is vertical and remains unchanging in po- 
larity (Fig. This results in having some concave and some 
convex lenses. 

In the straight sections, away from focusing fields, we add a 
weak bias field (about one kV/m for the drift spaces) to main- 
tain the quantization axisjl2t|. 



field. The different \M\ states (or M states as in Fig. Q) have 
different Stark shifts and hence experience different forces 
in an electric field gradient. If the molecules in a ring re- 
peatedly enter regions of weak and direction-changing elec- 
tric fields, transitions between the different M components 
will take place (Majorana transitions), leading to a loss of the 
molecules [pi] |l2| , p~3| ] . Our method of avoiding this problem 
is described below. 



C. Avoiding weak-field regions 

In the bend sections, which use electrodes with zero electric 
field near their center, the centripetal force keeps the trajectory 
of molecules in a strong electric field at the inner side of that 
zero (Fig. |]c). And, as did Crompvoets et al.[@], we take 
advantage of the combined horizontal and vertical focusing. 

Focusing lenses in the straight sections (Fig. ||a,b) are sex- 
tupolar with a dipole field added to avoid zero field. The mul- 
tipole coefficients of the electrodes, their gaps, and electric 
field strengths on orbit, are listed in Table II. With the added 
dipole field, these lenses focus in only one transverse direction 
while defocusing in the other and are used in an alternating- 



D. Dynamics calculation 

The force on a molecule in one direction is given by Eq. 
|l[ The basic formulas for motions of the molecules are de- 
rived analytically without linearization. Then the beam optics 
and dynamics are numerically calculated and optimized by us- 
ing a modified version of the program Goemon[|l5| that was 
originally developed for relativistic charged particles. The 
molecules were numerically tracked through the synchrotron 
for 400 turns (15 s) to determine the beam's dynamic stability. 

The numerical integration was done in the time domain be- 
cause the speed of a molecule varies as a function of the elec- 
tric field. The overall dimension of the lattice is chosen to 
balance easily attainable electric fields and high transverse ac- 
ceptance with compact (desk top) size. 



1. Potential energy 

The potential energy, W, in an electric field, of the \J,K >= 
1 1,1 > level of 14 N 2 H 3 (Fig. |) is given, in units of Joules 
(written out to avoid confusion with the quantum number "/") 



4 



by: 



W = ± 



Ci - C 3 E Z 



where Ci = 5.26 x l(T 25 Joules is half of the \J,K >= |1, 1 > 
inversion splitting (Fig. ||), C2 = 2.52 x 10~ 30 JoulesV _1 m, 
C 3 = 1.78 x 10" 38 JoulesV" 2 m 2 and £ is in V/m. The terms 



involving C\ and C2 are taken from Ref.[ 10 1 to which we add 
a second order term to account for mixing of the \J,K > = 
|2, 1 > state in strong electric fields [fiT 



4 

to 



■ 

- (a) "J 


ggg„ „ 








■■■■■■■■■■■ ■ 


lii" 


11 



■j 

(b) j 












III* 


11 



-4 4 
x (mm) 



-4 4 
x (mm) 



2. Electric field gradient 

Within an electrode, the electric field, as a function of dis- 
tances x and y from the reference orbit, is calculated from 
a scaler potential \\t by E = ~V\\t(x,y) and \\i is taken to be 
uniform in the longitudinal (z) direction within the electrode 
length. The fringe fields, at the ends of the electrodes, where 
the electric field is changing in the z direction of motion pro- 
duce additional transverse forces. These are evaluated sepa- 
rately, and are described in section II D 3. 

In a straight (focusing/defocusing) electrode, the potential 
is a combination of dipole and sextupole terms given by 
\j/ s = — Eo[y + Aj,(x 2 y — iy 3 )] where Eq is the central elec- 
tric field, A3 is the sextupole coefficient (see Table II) and the 
dipole coefficient has been set equal to 1 . 

In a bending electrode of constant radius p, \|/ is a combi- 
nation of dipole, quadrupole, and sextupole terms given by: 



-E b 



y + yB 2 \n 1 



-B 3 J [k(p+x)]smh(ky) 



Relating this to the Cartesian multipoles, we find: Et, = 
E [l +2A 3 //fc 2 ], B 2 = pA 2 [k 2 /k 2 + 2A 3 ], where A 2 is the 
quadrupole coefficient, 



#3 = - 



1 



2A 3 



kJ (kp) k 2 + 2A 3 ' 



Jo and J\ are the th and I st order Bessel functions and the 
value of pk is the first root of Ji(kp) = 0. An approximation 
for useful for comparison with \)f s up to third order is 



= -E 



y+A 2 xy+A 3 



1 - 



A 2 
2pA 3 



2 y 

x y 

y 3 



3. Fringing fields 

To evaluate the transverse vertical focusing forces in the 
fringing regions of the elements, a two - dimensional numer- 
ical calculation, using the geometries of the electrodes, was 
carried out to find the electric field as a function of z on the 
midplane. The magnitude of the fringe fields are shown in Fig. 
^. The focusing action was then calculated analytically from 
the derivatives of this field, evaluated as equivalent thin lenses 
at the ends of the focusing, bending, and bunching electrodes, 
and included in the calculation of trajectories around the ring. 



FIG. 7: Dynamic aperture for on-momentum molecules: the start- 
ing coordinates for the molecules that survive 400 turns (a) without 
gravity and (b) with gravity. Dynamic aperture is the area, in the 
transverse plane, at the center of the long straight section, occupied 
by the molecules. 



4. Hamiltonian 

All optics, orbit, and tracking calculations were carried out 
using a second-order symplectic integrator and the Hamilto- 
nian: H — Ho + W(x,y). In straight regions, Hq = j^(P 2 + 
P 2 + P 2 ) where m is the mass of the molecule, P x = mv x , 
P y = mVy, and P z = mv z is the momentum, with v x , v y , and 
v z the velocities in Cartesian coordinates. Conservation of 
energy requires that P z changes when the molecule passes 
through a field gradient at the ends of the electrodes. 

In a bending region, the Hamiltonian Hq becomes: 



1 



Hq = — 
2m 



(p+x) 2 



where Pq is an angular momentum, p is a reference bend 
radius, x is a radial excursion with respect to p, and P z = 
(Pe/p+x) 2 . 



III. SYNCHROTRON PERFORMANCE 
A. Dynamic aperture 

The survival of molecules, tracked through 400 turns, is cal- 
culated to determine the usable dynamic apertures, a x and a y . 
Fig. [7| shows a scatter plot of the starting coordinates x and 
y at center of the straight section of those molecules that sur- 
vive the 400 turns. From these we calculate the acceptances, 
e, as l/7t times the areas in displacement-angle phase space, 
e = a 2 /p. The dynamic apertures, values of beta at center of 
straight section, and acceptances are given in Table III. In Ta- 
ble III the circulation period is 38 ms whereas 3.357/90=37.3 
ms, reflecting the 3% to 4% reduction in velocity when the 
molecules enter the electric fields. 

The dynamic aperture of off-momentum molecules deter- 
mines the momentum acceptance. Fig. || shows the effect 
of momentum on the dynamic aperture, and indicates an ac- 
ceptance of about -3% to +1.7%. We shall see later that this 
reduces to about ±2% when the bunchers are operating. 
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FIG. 8: Dynamic aperture as a function of momentum variation from 
the nominal tune for x (solid points) and y (open points) with gravity 
but without bunching. The molecules in Fig. M appear here at dP/P = 
0. 



FIG. 9: Closed orbit displacement (COD) from the reference orbit 
due to gravity. The swings of the COD in the bending sections are 
similar to the curve of the horizontal dispersion because the vertical 
and horizontal betatron tune advances are both 4lt. 



TABLE III: Synchrotron operating parameters 



Circumference (m) 


3.357 


Circulation period (s) 


0.0380 


Velocity in free space (m/s) 


90.0 


Beta horizontal a : p A (m) 


1.264 


Beta vertical pVy (m) 


0.513 


Horizontal dispersion (m) 


0.001 


Horizontal tune: v x 


5.250 


Vertical tune: v v 


5.200 


Chromaticity - horizontal: ^ A 


-0.0885 


Chromaticity - vertical: £ y 


-0.0942 


Momentum compaction: a 


-0.99 


Dynamic aperture - horizontal: a x (mm) 


6.5 


Dynamic aperture - vertical: a y (mm) 


6.0 


Acceptance - horizontal: e A (mm - mr) 


35 


Acceptance - vertical: e v (mm - mr) 


71 


"At the center of the long straight section. 



The betatron tunes, as a function of momentum, are found 
by tracking the motions of a molecule over 512 turns. The 
chromaticities, are ^ = -0.0885 and Cy = -0.0942. The ef- 
fect of momentum-deviation upon the circulation period T, ex- 
pressed as the momentum-compaction factor a in Ar /7b = 
aAP/Po, is strong. The value of a is -0.991; this will result in 
prompt debunching if the bunching voltages are turned off. 



B. Effect of gravity 

The effect of gravity is visible in the orbit because the ve- 
locity of the beam is low. The vertical phase advance of 4n 
(same as the horizontal phase advance), is favorable in reduc- 
ing the effect of the gravity force on the vertical closed orbit. 
Fig. ^| shows the vertical closed-orbit displacement (COD) 
due to gravity. It ranges from -0.13 mm to 0.09 mm and is 
small enough to be neglected in the case of linear optics calcu- 
lations. In the case of dynamic aperture calculations we have 
confirmed that gravity changes only the details of the dynamic 
aperture at its edges and does not affect its core portion (Fig. 

0). 
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FIG. 10: Buncher timing for an on-energy molecule and a faster 
molecule. Changing the buncher timing may be used to change the 
circulating beam energy. 



C. Bunching and deceleration 

Bunching action is provided by six sets of parallel-plate 
electrodes, each 20 mm in length with 10 mm half-gap, placed 
in low dispersion regions around the ring (Fig.[j], H). The volt- 
age on the plates is pulsed in a triangular waveform that pro- 
duces a maximum electric field of 1.0 MV/m for on-energy 
molecules. Each buncher is synchronized to the on-energy 
molecules (Fig. |Io|) so that the energy lost upon entering the 
buncher is matched by the energy gained upon exit. A higher- 
energy molecule reaches the buncher at an earlier time when 
electric field gradient is higher upon entrance and lower upon 
exit, resulting in a net energy reduction for a molecule in a 
weak-field seeking state. 

The buncher frequency is 1.58 kHz and allows 60 bunches 
around the ring, spaced 56 mm apart. To prevent Majorana 
transitions the voltage is biased to avoid negative fields. The 
rate of change of the electric field in the bunchers can be made 
as large as about 7 GVm^'s -1 . The synchrotron tune is pro- 
portional to the square root of the rate of change and is 0.92 at 
our reference value of 7 GVm~ 1 s~ 1 . 

The lattice is optimized for a nominal velocity of 90 m/sec. 
After stacking at this velocity, the beam can be decelerated by 
synchronously changing the electrode and buncher settings. 
This provides an opportunity to scan the velocity during an 
experiment. The betatron tunes stay constant but the straight 
section becomes dispersive. At 63 m/sec, the horizontal trans- 
verse acceptance reduces from 35 to 13 mm-mrad and the mo- 
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FIG. 11: Scatter plot of an initially injected ±158/JS-long pulse of 
molecules with a ±2% momentum spread after 1 - 400 turns. The 
velocity spread of ±1.8 m/s represents an energy spread in the mov- 
ing frame of about ±0.4 K. 
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FIG. 12: Last elements of a linear decelerator showing a group of 
four decelerating electrodes and one focusing lens. The lengths of 
the decelerating elements and focusing lenses decrease to match the 
velocity of the decelerating molecules. 

seeking state in electrostatic traps. They find that the loss rates 
can be significant, in some cases precluding successful evap- 
orative cooling, and are influenced both by the electric fields 
and the size of the electric dipole moment. 

Thus, the use of a synchrotron storage ring with high den- 
sities of molecules in a weak-field seeking state will need to 
confront this issue. Obviously, the alternative would be to use 
molecules in the strong-field seeking rotational ground (/ = 0) 
state. A synchrotron storage ring for molecules in the J = 
state will be discussed in a future paper. 



IV. SYNCHROTRON STORAGE RING SYSTEM 

The overall performance of a molecular synchrotron stor- 
age ring also depends upon the beam delivered by the source, 
decelerator and injection line. To model performance and de- 
termine the stored beam from the complete system we also 
modeled a decelerator, injector and source. 



mentum acceptance from 2.0% to 1.2%. 



D. Collisional losses 

So far, we have only examined storage ring losses associ- 
ated with the limits of the dynamic aperture and momentum 
acceptance of the lattice. A real storage ring will also have 
losses due to elastic and inelastic scattering of the molecules. 
Scattering by room-temperature background gasses is the ma- 
jor elastic scattering contribution and of the gasses likely to be 
present in the ring, xenon and ammonia have the largest scat- 
tering coefficients. Assuming hard sphere binary collisions 
between the 14 N 2 H3 molecules and 300 K xenon (leaking in 
from the source), the mean time between collisions is 15s at a 
(xenon) pressure of 1 x 10~ 7 Pa (7.5 x 10~ 10 Torr), where we 
have used the equations and collision diameter in Ref . [p^| . 

For inelastic scattering between the 14 N 2 H3 molecules, 
their relative kinetic energy spread of about ± 0.4 K (±5.5 x 
10~ 24 J) is small enough to prevent any significant excitation 
to higher rotational states. A molecule can however collision- 
ally relax to a strong-field seeking state, causingit to be lost 
from the storage ring [^[|. Bohn [ |l8| ] and Kajita [ jl9| have cal- 
culated inelastic collision rates for molecules in a weak-field 



A. Linear decelerator 

Slow molecules, suitable for injection into a storage ring, 
can be produced by time-varying electric field gradient de- 
celeration, by mechanical cancellation of the molecular beam 
velocity [^o[, and possibly by buffer-gas cooling [|TJ (without 
magnetic trapping). Time-varying electric field gradient de- 
celeration is the easiest of these for us to model and to match 
to the storage ring. It has been used by Bethlem et al. [ p2[ [0| 
and Crompvoets et al. [0] to decelerate 14 N 2 H3. 

Our model linear decelerator takes the 310 m/s (115 K ki- 
netic energy) output of a 14 N 2 H3 - seeded xenon pulsed-jet 
source (room temperature reservoir) and decelerates it to 90 
m/s (9.75 K kinetic energy). The number of electrodes is set 
by the decrease in kinetic energy in each electrode (equal to 
the change in potential energy of the molecule entering the 
electric field). We use 79 decelerating electrodes, decreasing 
in length from 48 mm to 8.75 mm effective length to keep 
the transit time of the molecules constant through each elec- 
trode. We choose the overall length of the decelerator, 3.4 m, 
comparable to the 3.36 m circumference of the storage ring, 
to balance high velocity acceptance with compact size. 

The electric field is nearly a square wave 158 /js long with 
a repetition rate of 3. 16 kHz and the maximum electric field is 
8 MV/m in a 1.75 mm half-gap and is the same for all decel- 
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TABLE IV: Injection matching fields 



-4 4 
displacement (mm) 



-4 4 
displacement (mm) 



FIG. 13: Boundaries of injected-beam emittances (solid lines) and 
storage ring acceptances (broken lines) in the (a) horizontal and (b) 
vertical planes. 



erating electrodes. After the molecules enter the electric field 
it drops to nearly zero so that the molecules that exit the elec- 
trode must relinquish kinetic energy to enter the field in the 
next set of electrodes. This is done either by having succes- 
sive electrodes 180 degrees out of phase or by using only ev- 
ery other bunch. The electric field does not return completely 
to zero and the horizontal focusing elements (see below) have 
their electric field in the same direction as the decelerating 
electrodes to minimize Majorana transitions. 

Bunching, as the molecules decelerate, is accomplished by 
having the electric field (in each bunching electrode) decrease, 
linearly over the 158 /js from 8 MV/m to 7.76 MV/m so that 
the fastest molecules, arriving early, receive the most decel- 
eration. This results in a ±1% momentum spread at 90 m/s 
which falls within the momentum acceptance of the storage 
ring. 

The spacing between individual decelerating electrodes re- 
mains constant. This fixes the fringe field which provides 
vertical focusing. The decelerating electrodes are grouped 
in sets of four as shown in Fig. [12|. After each quadruplet 
of decelerating electrodes is placed a horizontal focusing ele- 
ment with length appropriate to keep the molecules in phase 
with the next decelerating quadruplet. (One section has only 
three decelerating electrodes.) The 20 focusing elements all 
have the same field but their focusing strength increases as 
the molecules slow. The lattice is, in the vertical direction, 
FFFFD and in the horizontal direction OOOOF. The overall 
length of the decelerator is 3.4 m. 

Upon exit from the decelerator, bunches are 7 mm long with 
a vertical and horizontal half-width of 1.30 mm and 1.61 mm, 
respectively. The vertical and horizontal betas are: P v = 56.3 
mm and p/, = 86.7 mm corresponding to vertical and hori- 
zontal emittances both of 30 mm-mr. Being less than the 7 1 
mm-mr and 35 mm-mr vertical and horizontal acceptances of 
the synchrotron storage ring, this sets the overall transverse 
acceptance of the entire system. 



B. Injector 

A bunched beam of molecules is injected onto the closed or- 
bit of the synchrotron at the downstream end of a straight sec- 
tion. Starting at the exit of the decelerator, this beam is guided 
along a trajectory by electric fields that focus transversely and 



Electric field (MV/m) 
Dipole coefficient (m 
Radius (m) 
Arc length (rad) 
Phase advance (rad) 
Focus parameter, P (m) 



Sector 1 

3.37 

158.5 

-0.20 

0.397 

1.57 

0.05 



Sector 2 

4.0 

38.2 

0.69 

0.426 

1.55 

0.19 



transform its size and divergence to match the vertical and 
horizontal acceptances of the storage ring. 

The last element must be a pulsed deflecting electrode that 
turns off after the bunch or bunches have entered. Several 
different injection protocols are available: a string of 158 /js 
pulses can be injected into the ring in less than the circulation 
period of 38 ms, or single bunches may be filled at any time by 
switching the pulsed deflecting electrode during the interval 
between bunches. 

The injector guide field consists of two bending sectors, 
each an arc of about 7t/8 rad, arranged as shown in Figjj]. In 
each bending sector the field is configured to provide equal 
horizontal and vertical focusing. The phase advance of trans- 
verse motion is near n/2 in a sector; thus, it exchanges an- 
gle for displacement and changes the ratio of angle to dis- 
placement. Focusing strengths are adjusted to provide the 
required match between injecting beam and ring acceptance 
The boundaries of the beam phase space are shown in Fig. 13 

It was possible to choose the first sector to have radius and 
strength equal to that of the bend in the storage ring. The sec- 
ond sector must have larger aperture and radius and it must 
be pulsed to zero at the end of injection. Parameters of the 
injector guide field are shown in Table IV. The sum of the 
two inverse curves was made near zero so that the line of the 
decelerator is about parallel to the straight section of the syn- 
chrotron storage ring. 



C. Source and intensity 

To estimate the number of molecules that can be deceler- 
ated, injected, and stored, we assume 14 N 2 H 3 source condi- 
tions similar to those reported by Crompvoets et al. [|7| and by 
Bethlem et al.[^]: a pulsed jet source of 0.8% 14 N 2 H 3 seeded 
in 152 kPa (1140 Torr) of xenon (reservoir temperature 300 
K) exiting through a circular 0.80 mm dia. orifice into vac- 
uum. Following Crompvoets et al. [^J, we assume that 15% of 
the molecules entering the decelerator are in the desired state. 
From the formulas in Miller|Q, we find a flux of 4.4 x 10 18 
14 N 2 H3 molecules sr^'s -1 in the desired M state. 

The decelerator's transverse emittance of ±30 mm-mr and 
momentum spread of ± 1 % at 90 m/s set a transverse accep- 
tance from the 310 m/s jet source beam of ±9 mm-mr and mo- 
mentum spread of ±0.29%. If we assume an initial Gaussian 
velocity distribution with a mean of 310 m/s and a standard 
deviation of ±20%, approximately 1.1 % of the molecules 
from the jet source fall within our decelerator momentum ac- 
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ceptance. With a pulse length of 158 /js, the intensity of the 
decelerated beam is roughly 6x10 molecules/pulse. 

The horizontal acceptance of the storage ring matches the 
emittance of the decelerator and the vertical acceptance and 
momentum acceptance of the storage ring is about twice the 
emittance of the injected beam, so all 6 x 10 s molecules/pulse 
should be captured and stored in a single bunch. If all sixty 
bunches are filled (over multiple turns) the total stored beam 
is 3.6 x 10 10 molecules and the circulating flux is 9.5 x 10 11 
molecules/s. 
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